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M atrix programming of 

differential equations for electronic 
analog computers organizes the equations, 
thereby minimizing chances for errors; 
simplifies the scaling of coefficients and 
variables to suit the characteristics of the 
available machine; and in many instances 
permits the reduction of the number of 
“sign changers.” The matrix method 
furthermore helps safeguard against the 
appearance of extraneous terms in the 
solutions of certain systems of differential 
equations. These very practical attri¬ 
butes of the matrix method of program¬ 
ming save computing time and equipment. 

The key concept in the theory of matrix 
methods of programming electronic ana¬ 
log computers lies in the establishment of 
the 1-to-l reciprocal correspondence be¬ 
tween the matrically formulated differ¬ 
ential equations and the synthesized 
electronic computer networks. 1-3 The 
generality of the matrix method leads not 
merely to one or two, but to a large group 
of electronic synthesis networks all yield¬ 
ing the solution to a given differential 
equation problem. 

The object of this paper is to indicate 


by specific examples arising in practice 
the advantages of the matrix program¬ 
ming of electronic analog computers. 

In the first example a wind-tunnel 
problem illustrates the organization, 
scaling of coefficients, and reduction of 
negative signs, thus minimizing the need 
for sign-changing amplifiers in a high- 
order system of differential equations. 

This is followed by an elementary ex¬ 
ample exhibiting the care which must be 
taken in programming certain systems of 
differential equations to preclude the 
appearance of extraneous roots in the 
determinantal equation, and thereby 
extraneous solutions to the equations. 
This particularly insidious phenomenon 
may arise from devotion to the common 
“solve for the highest derivative” routine 
for synthesis of computer networks. 
These extraneous solutions are some¬ 
times manifest as “drift” rather than in 
their true nature. Furthermore, the 
identical computer network sometimes 
yields “solutions” which are correct, and 
sometimes incorrect, depending upon the 
initial conditions. This is illuminated by 
a specific example. 


Elementary Matrix Transformations 

The objectives of some linear trans¬ 
formations of systems of linear algebraic 
equations and first-order differential 
equation Systems of particular usefulness 
in analog computers are: 1. changing 
the magnitudes or signs of the coefficients 
of the terms in the equations, and 2. 
scaling in magnitude or changing the 
signs of the prescribed functions or de¬ 
pendent variables. Of course, such 
changes as are made must yield solutions 
unchanged or easily related to the original 
problem equations. The need fojc these 
transformations is often dictated by the 
purely electrical characteristics of the 
computer .or recording equipment, such 
as the allowable output voltage range of 
the amplifiers. These elementary trans¬ 
formations are illustrated by the following 
example. 

Example 1. The design of a wind- 
tunnel control system using analog com¬ 
puter techniques requires synthesis of the 
equations describing tunnel performance 
and the experimental determination of 
desirable controller characteristics. The 
particular wind-tunnel considered in this 
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example may be adequately represented 
by the following differential equation sys¬ 
tem, in which P represents the pressures 
developed as a function of the prescribed 
variables, A and A’ 0 . (The differential 
coefficient d/dt is represented by d.) 

(1.90 +rf) P, - 0.423P< -1.59 P ac = 

-0.303^i t (1) 

- 0.433P,+(4.51 +<i) P 4 -1.87 P md - 

0.331P n -4.K>P p =23.7A,<,-1.29N, (2) 

- 58.1P,+(74.9+d)P«-30.8P P = 

155A r 0 (3) 

-1.01P 4 +(0.124+d)P md -7.52P, = 

0.302AT* (4) 

(0.309 +<f) P m i —0.159P m = 37.8AT 0 (5) 

-0.483P m< +(0.617+d)P m = 

—106..4 so —108..4 is (6)‘ 

-0.00514P,-0.091P«+(2.97+d)P p - 

0.973P, =0.173.^,,-2.18AT, (7) 

-15.6P md +(16.3+<f)P,=0 (8) 

151.P„-54.3P,+10‘P„/P p =0 (9) 

Valve openings which determine pressure 
are represented by A and are assumed to 
be controlled by functions of the differ¬ 
ences between desired pressures Pm») Ptt, 
or pressure ratio ( P n /P v ) 9t and respective 
measured values Pm) P*) ( Pn/Pp)' Con¬ 
troller characteristics may be represented 
functionally by the following equations 

^3» —/i(Fm*— Pin) (f0) 

•^30 = ft(Ps8 ~Pt) ( 11 ) 

Aa=M{Pn/Pp)s-(Pn/Pp)} (12) 

To illustrate application of the matrix 
technique to this problem, the synthesis 
procedure for this 8th-order system is 
described. Inclusion of controller char¬ 
acteristics involves an addition to the 
basic matrix equation and is not essential 
for the purpose of the present example. 

An equivalent matrix representation of 
the system of equations 1 through 9 is 
indicated in equation 13 (see below). 

The problem is to adjust the coefficients 
to reasonable values for use with a par¬ 
ticular computer. Since the dependent 
variables are not known a priori, the term 
reasonable initially involves a “best 
guess” that tends to keep both the mag¬ 


nitude of the coefficients and dependent 
variables within acceptable limits. A 
convenient method is to apply transforma¬ 
tions that attempt to establish coefficients 
approaching unity. For the majority of 
electronic analog computers, the step 
between the final mathematical form and 
the synthesis then involves the introduc¬ 
tion of a factor of 10 “ 6 in all terms, thus 
permitting immediate interpretation in 
terms of computer component values 
expressed in microfarads and micromhos, 
and computer variables expressed in volts 
and microamperes. 

In equation 13 premultipliers affecting 
rows 3, 8, and 9, as follows: 

row 3—divide by 100 
row 8—divide by 10 
row 9—divide by 100 

result in reasonable values for all terms 
except that in row 9, column 9, which re¬ 
mains beyond normal limits. For cases 
of this type application of postmultipliers 
affecting columns in the coefficient matrix 
and also the dependent variables achieves 
the desired result. In particular, applying 
the column transformation 

column 9—divide by 1,000 

completes the coefficient adjustment 
necessary in equation 13. 

Treatment of coefficients having an 
associated negative sign conventionally 
requires an additional amplifier to act as 
a sign changer. Reduction of the number 
of such sign changer is of advantage from 
both equipment and stability standpoints. 
Using the matrix formulation, transforma¬ 
tions changing sign are easily applied and 
the effects of each transformation relative 
to the reduction of minus signs is imme¬ 
diately evident. Although unnecessary 
mathematically, electronic network sta¬ 
bility usually dictates that transforma¬ 
tion of rows and columns are made in pairs 
thus maintaining principal diagonal terms 
with positive sign. Desirable transforma¬ 
tions for equation 12 change the signs 
of: 

row 1, column 1 
row 4, column 4 
row 6, column 6 
row 7, column 7 
row 9, column 9 
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-1.59 

0 

0 
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— 0.303.4 u 

-0.433 

(4.51-H0 
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-1.87 
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-0.0331 

-4.16 

0 
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pi 


23.7.4 »o-l. 29 A r 0 

-58.1 

0 

(74.9+d) 

0 

0 

0 

-30.8 

0 

0 

Pac 


155.A r 0 

0 

-1.01 

0 

(0.124+J) 

0 

0 

0 

-7.52 

0 

P md 


0.302M, 

0 

0 

0 

0 

(0.309 -M) 

-0.159 

0 

0 

0 

P mi 

= 

37.8A T o 

0 

0 

0 

0 

-0.483 

(0.617+d) 

0 

0 

0 

Pm 


— I00..4jo-X08..4ji 

-0.00514 

0 

-0.091 

0 

0 

0 

(2.97 -M) 

-0.973 

0 

Pv 


0.173.4, 2 -2.18.V 0 

0 

0 

0 

-15.6 

0 

0 

0 

(16.3+d) 

0 

Ps 


0 

0 

0 

0 

0 

0 

0 

151. 

-54.3 

10 5 

L Pn/Pp - 


0 


(13) 
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The number of sign changers has now 
been reduced from 17 to one. 

Results of the application of both scale- 
factor and sign-change transformations 
are indicated in equation 14 (see pre¬ 
vious page), which is equivalent to equa¬ 
tion 13 but more applicable to analog 
computer synthesis. 

If equation 13 is represented as 


[M*o] 

= bol 

(15) 

and equation 14 as 



1 = Oil 

(16) 


Extraneous Solutions 

The elementary solve-for-the-highest- 
derivative technique of computer pro¬ 
gramming, although restricted in form, 
usually gives correct solutions, partic¬ 
ularly for single wth-order differential 
equations. In systems of differential 
equations, however, great care must be 
exercised to avoid extraneous solutions. 

It is knowm that differentiation of equa¬ 
tions in a system may introduce additional 
roots in the determinantal equation, and 
thus alien solutions. This constraint is 


equations 17 through 19 represent 
elementary transformations. 
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-1 0 
0 1 
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0 
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0 

0 

0 

-1 
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0 

0 

0 

0 

0 

0 

0 

0 
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(17) 


(18) 


(19) 


In view of their elementary character, the 
matrix premultipliers and postmultipliers 
in equations 17 through 19 need not be 
written down explicitly, although the in¬ 
formation they contain, especially in the 
last, must be recorded. 

The electronic network equation corre¬ 
sponding to equation 14 is indicated in 
equation 20, in which the coefficients now 
represent micromhos and microfarads, 
and the resultant ideograph is shown in 
Fig. 1. 


gram a computer according to the cur¬ 
rently accepted methods which require 
that each equation be written in such a 
manner that a derivative of a dependent 
variable appears to be a known function. 
Then, assuming the dependent variables 
are known, successive integrations are 
performed, functions combined, and 
eventually “fed back,” completing the 
synthesis network. This common tech¬ 


nique is not applicable to many systems of 
equations and the differentiation required 
in attempting to apply the method may 
easily yield erroneous results. Example 2 
is an elementary illustration of this type 
of programming error. 

Example 2. The following system of 
differential equations 

d 2 x +d 2 y — d 2 z +dy +.r +y = 0 

dy+z = 0 (21) 

2dx+dz+z = 0 


has the determinantal equation 

X(X+1)(X 2 +1)=0 (22) 

and the general solution 

• (C3~ C 2 ) , (Cj+C|) . 

x =-cos t -sin l — c K 

2 2 

y = C\e~ 1 — ci sin t +c 3 cos t +c 4 (23) 

z = C\e~ t +c 2 cos /+c 3 sin t, 


in which the four constants specified in 
equation 23 are determined by the con¬ 
ditions 


x(0) = 


to — c 2 ) 

2 


— C 4 


y( 0) =Ci+C3+C4 

z(0)=d+c 2 (24) 


dx( 0) = - 


to+£s) 

2 


It may appear desirable to differentiate 
the second and third equations of equa¬ 
tion 21 yielding a set of equations which 
may then be programmed without diffi¬ 
culty using the conventional procedure. 
The new system reads 

d 2 x -\-d 2 y—d 2 z +dy +x -f -y = 0 

d 2 y+dz = 0 (25) 

2d 2 x -\-d 2 z +dz — 0 

and has the determinantal equation 
X 3 (X+1)(X 2 + 1)=0 (26) 

and the general solution 

to — £ 2 ) , toTc 3 ) • , 

x — -cos t --— sin t — 

2 2 

to+c«) ~ c & 

y = Cie~ l —C 2 sin t+c 3 cos /+c 4 +c 6 / (27) 

z = Cie~ l -\-C 2 cos t+c 3 sin t+c 5 
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Fig. 1. Ideograph of equation 20 


The six constants in equation 27 are 
related to the initial conditions by the 
relations 

/r\\ ( C 3-C 2 ) 

*(0) = —o- C *-Cs 


y(0) — Ci+C3+C4 


z(0) =ci+c 2 +r5 


dx(0)= — 


fe+fa) 

2 


— Ce 


dy( 0) = -ci-c 2 +ce 
dz(0) = —CiTcj 


(28) 


dicates the appearance of extraneous 
constant and linear terms in the latter. 

In the general situation, the solutions 
to the differential equation system are not 
known explicitly (else why use the com¬ 
puter?). This leads to the following 
interesting situation. Assume that the 
conventional technique has been followed 
and the computer set up according to the 
differentiated system equation 25. Now 
suppose the transient response to the 
initial conditions 

*( 0 ) =20 

y(0) =z(0) =dx(0) =0 (29) 


*=-(20 cost +40) 

O 

y =-(20e~ 1 +20 cos *+20 sin t— 40) (30) 

o 

z = -(20e~ t —20 cos /+20 sin /), 
o 

i.e., either system yields the identical 
solution. 

Next, consider the solution for the 
initial conditions 

x(0) =y(0) =z(0) =0 

dx( 0) =20 (31) 


Comparison of the solution in equation 
23 with that in equation 27 clearly in¬ 


is desired. Then the solutions given by The solution in equation 23 to the original 
equations 23 and 27 both read problem of equation 21 now yields 
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Fig. 2. Ideograph of equation 34 



Fig. 3. Ideograph of equation 38 


x = 20 sin t 

y = 20e~*+20 sin /-20 cos t (32) 

z=20e~ 1 —20 sin /—20 cos t 

whereas equation 27 yields the distinctly 
different solution 

x =*-(20 sin /+40 t) 

o 

y =~( —20e“ r -f20 sin /—20 cos t+ 

o 

40-40/) (33) 

z = -( — 20e~'—20 sin t—20 cos /+40) 

O 

Recalhng that the computer has been 
set up according to the common tech¬ 
nique, i.e., according to equation 25, 
it is clear that the computing machine 
(which merely exhibits graphical plots) will 
in the first instance, equations 29 and 30, 
yield correct solutions to the ostensible 
problem in equation 21. In the second 
instance, equation 31, the computer 
exhibits solutions which do not satisfy 
the original problem in equation 21. 

Of course, many other combinations of 
initial conditions are possible; but the 
machine never divulges which are indeed 
the correct solutions to the original 
problem, and which are not. This must 
be ascertained by other means. In this 
instance, the unwary may easily be misled, 
for as shown, sometimes the correct 
solutions are obtained, and sometimes 
they are not. 

The matrix formulation of the problem, 
however, places in a clearer light the well- 
known mathematical requirements on 
transformations such that any synthesized 
computer network, if stable, yields the 
solutions to the original problem. This is 
indicated in the next section, together 
with oscillograms of computer solutions 
to the differential equation problem of 
equations 21 and 25 programmed by the 
matrix method. 

Matrix Programming 

The matrix programming of the differ¬ 
ential system of equation 21 has been in¬ 
dicated and explained in reference 1, and 
one version leads to the equation 
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which is illustrated by the ideograph of 
Fig. 2. It therefore suffices to synthesize 
the differentiated system of equation 25. 
According to the method outlined in the 
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the first reference, equation 25 may be written as 
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in which Xi, x 2 , and x 3 represent x, y, and z 
respectively. 

One first-order system equivalent to 
equation 35 reads 
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(36) 

Since this is likely to be unstable on exist¬ 
ing analog computers, see reference 2, 
equation 36 may be premultiplied by 
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resulting in the matrix equation 38 which 
may be synthesized without difficulty 
and is stable. 
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(38) 


Fig. 3 is an ideograph of the electronic 
computer synthesis network correspond¬ 
ing to equation 38, and may be compared 
with Fig. 2. 

Comparison of solutions exhibited by 
an analog computer for the two sets of 
differential equations corroborate the 
assertions of the previous section. 

Figs. 4(A) and (B) respectively show 
the computer solutions to equations 21 
and 25 programmed according to the 
matrix equations 34 and 38 for initial 
conditions of equation 29. As stated, the 
machine solutions coincide in this case. 

For the initial conditions of equation 
31, the oscillograms Figs. 5(A) and (B) 
correspond respectively to equations 32 
and 33. As predicted, there is now no 
correlation between the respective solu¬ 
tions exhibited by the computer. 

Resume 

The practical merits of the matrix pro¬ 
gramming of analog computers have been 
briefly demonstrated in this paper. The 
most elementary of matrix transforma¬ 
tions are an invaluable time- and equip¬ 


ment-saving process, as the first example 
indicates. 

With regard to the second example 
and the introduction of extraneous roots 
in a system of differential equations, so 
strikingly demonstrated by the oscillo¬ 
grams, again matrix methods are ex¬ 
tremely valuable. 


It is well known that to obtain equiv¬ 
alent differential equation systems, the 
determinant of the transformation matrix 
must be a constant. In particular, the 
system of equation 38 may be obtained 
directly from equation 34, if the latter is 
premultiplied by the matrix 
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The determinant of equation 39 is 3d*. 



(A) (B) 


Fig. 4. Oscillograms for initial conditions of equation 31 

A—Programmed according to equation 21 
B—Programmed according to equation 25 



(A) (B) 


Fig. 5. Oscillograms for initial conditions of equation 31 

A—Programmed according to equation 21 
B—Programmed according to equation 25 
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Clearly this is not a constant and accounts 
for the extraneous X 2 term in the deter- 
minantal equation 26, and additional 
terms in the solution equation 33. 

In many practical problems which con¬ 
tain equations that have been differ¬ 
entiated for programming purposes, the 
manner of excitation of the network and 
initial conditions are such that the terms 
corresponding to the extraneous roots 
should have zero coefficients. However, 
whether due to initial or residual charges 
on the capacitors, grid current, internal 
parasitic coupling, or inexact matching 
of components, a small constant or linear 


Discussion 

Louis A. Pipes (University of California, 
Los Angeles, Calif.): Professor Honnell is 
one of the most outstanding leaders in 
stressing the importance of matrices and 
matrix calculus in the solution of problems 
in engineering and physics. Through his 
efforts and those of others it is now well- 
known that the theories of elastic structures, 
electric circuits, wave propagation, and 
mechanical vibration may be formulated 
concisely by the use of matrices, and that 
practical results may be obtained by means 
of the theorems of matrix algebra. 

The development of modern instruments 
of computation has created new matrix 
problems. The utility of matrices and 
matrix algebra in the formulation and solu¬ 
tion of many physical problems by the use of 
digital computing machines is now quite 
well known and is standard practice. Pro¬ 
fessor Honnell and Dr. Horn are to be con¬ 
gratulated for demonstrating the power and 
elegance of matrix methods in systematizing 
the programming of physical problems in 
analog computing machines. 

It is to be hoped that when the simplicity 
and conciseness of the method piesented by 
the authors are fully recognized, it will 
become a standard procedure in the pro¬ 
gramming of physical problems to be solved 
by analog computers. 


Vincent C. Rideout (University of Wiscon¬ 
sin, Madison, Wis.): As more complex 
sets of equations are attempted on analog 
computers, systematic methods, both for 
scaling and for manipulation prior to setting 
up block diagrams, become more necessary. 
The matrix methods introduced in this 
paper provide a valuable guide to those con¬ 
fronted by such problems. Although pre¬ 
sented for equations representing linear 
systems, they are also of help in nonlinear 
system study, if linearized forms of the 
system equations are first examined. 

The introduction of extraneous roots by 
differentiation is clearly pointed out in the 
second example of the paper. However, it 
would be interesting to know if the equa¬ 
tions of this example could come directly 
from a physical system. Also, this set may 
be shown to require only four arbitrary 
constants by the fact that it can be reduced 
to a third-order equation in y with the intro- 


term may be excited in addition to the 
correct solution response. The presence 
of such linear terms in the exhibited solu¬ 
tion may then be interpreted as amplifier 
“drift,” rather than actually as an in¬ 
correct synthesis network. This is the 
situation which may easily arise as in¬ 
dicated by the second example. 

Finally, it should be clear that the 
application of matrix methods to analog 
computer programming are of far wider 
scope and generality than can be indicated 
in this brief introduction. But the prac¬ 
ticality of the matrix method is clearly 
established. If not employed, the most 

- ♦- 

duction of only one arbitrary constant in 
the process of this reduction. 


Henry M. Paynter (American Center for 
Analog Computing, Boston, Mass.)’ and 
Daniel H. Sheingold (George A. Philbrick 
Researches, Inc., Boston, Mass.): The 
authors’ use of the matrix approach is to be 
commended for the ease of scaling, the 
economy of amplifiers, and the inherent 
solution stability, all of which are direct 
consequences of a highly organized and well- 
thought-out formulation technique. The 
comprehensive treatment of the subject of 
redundant integrations and ambiguous solu¬ 
tions is very timely and should be of para¬ 
mount importance to all those seriously 
interested in the machine solutions of sys¬ 
tems of differential equations. 

It is also particularly interesting to note 
the compact matrix form of the amplifier 
diagram; one novel feature of the symbolism 
is the use of conductance measures rather 
than resistances which results in a direct 
relationship between coefficients in the 
equations and physical parameters in the 
machine. However, we might suggest that 
a logical advantage could result from placing 
the feedback conductor in such ideographs 
in a position analogous to the capacitor 
whenever an amplifier is used merely as a 
summer or inverter. It is recognized, 
nevertheless, that there are perhaps com¬ 
parable advantages to having all capacitors 
north west/southeast and all conductors 
northeast/southwest. 

The principal purpose of this discussion 
is to show that certain commercially avail¬ 
able computers permit direct translation 



Fig. 6. Front view of a Universal Linear 
Operator showing decade parameter settings 
and external connections 


economic operational use of the machine 
is not realized. This is the antithesis of 
the raison d'etre for the computing ma¬ 
chine in the first place! 
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from the original matrix form (and, of 
course, from the original equations) to 
machine interconnections and parametric 
settings. Such machines are composed of 
modular elements which are called, for con¬ 
venience, Universal Linear Operators. Fig. 
6 shows the front panel of one modular 
element, and Fig. 7 shows a bank of 12 such 
elements. The basic equation, which com¬ 
pletely describes the performance of each 
element is 

4 

e = £o + 10 m p ^ djCj (40) 

i = i 

where 

m = 0, 1, 2, 3 
n =0, 1 

l/p = /o‘( )dt (41) 

The signal e Q is adjustable to any positive or 
negative constant voltage within the output 
range of the element in steps of 0.1 volt; 



Fig. 7. A bank of 12 Universal Linear 
Operators 
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K5_l_ SET-UP DIAGRAM 
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Matrix Programming of Electronic Analo.i Comoutors 


Set up for equation 2Q 



SWITCHES NOT MARKED. SET AT O OR NEUTRAL 


D - Direct Sura 
I - Integrate 


Fig. 8 (left). 
Tabular setup for 
equation 20 


Fig. 9 (right). 
Tabular setup for 
equation 34 
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SWITCHES NOT MARKED. SET AT O OR NEUTRAL 


D - Direct Sum 
1 - Integrate 
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SWITCHES NOT MARKED. SET AT O OR NEUTRAL. 


D - Direct Sum 
I - Integrate 

and is adjustable to any positive or nega¬ 
tive value from 0 to 11.10 in steps of 0.01 by 
means of 3-digit decade settings. It is seen 
immediately that this operator may func¬ 
tion either as an integrator (n = 1) or as a 



equation 34 for 
Xi, x 2/ and x 3 


summer (n= 0), and that the m parameter 
permits, in effect, a choice among four 
values of feedback resistance or capacitance. 
The voltage e 0 then functions either as a 
directly additive constant in the case of 


simple summation, or as an initial condition 
in the case of integration. 

The only fundamental distinction affect¬ 
ing problem formulation between the 
operators described in the foregoing and 
operational summing points in the authors’ 
matrix diagrams is in the location of the 
inverting amplifiers. The operational ele¬ 
ments described here have outputs which are 
of one sense only, i.e., ’‘direct,” or positive. 
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The choice of signs is made via a switch at 
each input. This somewhat simplifies 
programming formulation because all vari¬ 
ables may be defined as positive and be 
given sensed coefficients where operated 
upon. Although more than the minimum 
number of amplifiers might be used for a 
given problem configuration, the fact that 
they are inherently available in each oper¬ 
ator means that there is, in reality, no 
“waste” of amplifiers. 

We show in Figs. 8 through 10 tabular 
forms which completely define the computor 
problem setup including: 

1. interconnections 

2. parameters 

3. initial conditions 

and furthermore do so without reference to 
internal electronic functioning. Thus, in¬ 
dividuals without prior knowledge of the 
nature of analog computing circuitry can 
successfully set up such programs directly 
from the matrix equations, and entirely non¬ 
technical personnel can set up and run the 
machine directly from the tabular forms 
depicted. 

Fig. 8 gives the tabular setup for equation 
20 and Fig. 1 (for the homogeneous case). 
Note that the module numbers correspnd 
one to one with the problem variables (the 
columns in the matrix). Note also that 
when more than four inputs are required as in 
the second row of equation 20, additional 
modules can be employed; in this case only 
one extra module (number 10) was required. 

Figs. 9 and 10 correspond to equations 34 
and 38 and Figs. 2 and 3 respectively. 
Note again that the parameters can be set 
directly from the matrix itself. Fig. 11 
gives solutions of the second example under 
the initial conditions of equation 29, com¬ 
parable to those shown in Fig. 4 for x it x^, 
and x%. If these three variables and their 
rates of change must be displayed then six 
operators must be used. However, if only 
the variables themselves are required, the 
minimum number of four operators may 
be used. 

It may be of interest to note that the 
actual setup time requred by an individual 
reasonably experienced in the use of the 
equipment amounted to 4 minutes for equa¬ 
tion 20, and 2 minutes apiece for the other 
two examples. The actual running time for 
each case being under 1 minute, the total 
time required for the solution of all three 
equations was less than 15 minutes. 


V. H. Disney (Armour Research Founda¬ 
tion, Chicago, Ill.): The authors have de¬ 
veloped a very clear and useful system for 
the orderly manipulation of the equations 
to be mechanized on an analog computer. 
Description of the equations in matrix form 
permits an easily checked layout of the 
programming connections which should do 
much to reduce the possibility of error. It 
should be noted, however, that only the 
scales of coefficients and dependent variables 
are changed by the procedure described. It 
is suggested that this system be extended to 
include change of time scale as well. If this 
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were done the question arises as to what 
order of scale change in time, coefficient, and 
variable should be made for optimum pro¬ 
gramming. Guides might be developed for 
selecting reasonable values of the scale 
changes. 

The authors have also described the diffi¬ 
culties that can arise when equations are 
differentiated. This is not a widely used 
technique. I should like to ask the authors 
if they would recommend avoiding this 
practice altogether, or whether there are 
means to test the solutions for validity. 


R. E. Horn and P. M. Honnell: We are 
keenly appreciative of the favorable view 
of our work expressed by Dr. Pipes, for he 
is among the handful of mathematically in¬ 
clined scientists in the entire world who have 
pioneered and steadfastly pursued researches 
for almost a quarter of a century in the ap¬ 
plication of matrix methods to engineering 
problems. This was long before the elec¬ 
trical engineering profession generally had 
any cognizance of the existence of these 
mathematic methods, much less of their 
application to electical engineering prob¬ 
lems. 

The writings of Dr. Pipes have, over the 
years been a source of instruction, guidance, 
and motivation to all of us who now pursue 
matric methods in our own specialties of the 
moment. In particular, we hope to be able 
to apply the results of the researches of Dr. 
Pipes, in the difficult area of variable and 
nonlinear systems, to the electronic com¬ 
puting techniques we describe; needless to 


With regard to Dr. Rideout’s question 
concerning the origin of the equation in the 
second example, the basic form of the equa¬ 
tion was obtained directly from a system of 
equations representing an aircraft hydraulic 
servo system. In the original simulation of 
this problem, the extraneous linear term was 
introduced inadvertently using conventional 
techniques and caused considerable diffi¬ 
culty in the manner that is described in the 
paper. 

With reference to the discussion by Dr. 
Paynter and Mr. Sbeingold the following 
applies. One of the objectives of the 
authors’ paper was to indicate the applica¬ 
tion of the matrix approach to existing com¬ 
puters. By presenting such a complete dis¬ 
cussion of the application to a particular 
commercial machine, the discussors have 
placed additional emphasis on this aspect 
of the use of matrix methods. 

Because of their generality, the matrix 
programming methods should also aid con¬ 
siderably in obtaining the fullest utilization 
of other existing computers. 

We are glad to have Mr. Disney bring out 
the question of change-of-time scaling. 
This is again most concisely stated in terms 
of premultipliers and postmultipliers. Con¬ 
sider, for example, the change in time scale 
represented by 

r/k—t and/(r/&) =/(/), 
where 

/ = original independent variable 
r = transformed independent variable 
k = constant of proportionality 

as applied to equation 22 of reference 3, 
namely 
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With appropriate premultipliers and post¬ 
multipliers (omitted herein), this yields the 
desired expression 
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wherein 


say, these are matric methods! It must 
give Dr. Pipes great satisfactoin to see his 
prescience so abundantly substantiated. 

The comments of Dr. Rideout concerning 
the application of these methods to systems 
of nonlinear equations are particularly per¬ 
tinent. It is in these areas that a large 
percentage of problems of current interest 
lie—complex systems of equations in Which 
the systematic matrix approach is partic¬ 
ularly helpful. The authors plan to con¬ 
sider these aspects in more detail in their 
future publications. 


Dr n =lXdr n , 

d T n is written for d n /dT n , and, of course, D = 
(1 /k)D r . 

With respect to Mr. Disney’s last point, 
we defer to the mathematicians, who teach 
that equivalent differential systems can be 
obtained only if the determinant of the 
transformation matrices are constant. This 
permits the appearance of the derivative 
operator among the elements of transforma¬ 
tion matrices, subject to the aforementioned 
restriction. It does rule out the differentia¬ 
tion of equations per se. 
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